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ORBIT MODELLING 


FOR SATELLITES USING THE 
NASA PREDICTION BULLETINS 

N. L. Bonavito 
D. W. Koch 
G. A. Maslyar 
J. C. Foreman 

ABSTRACT 

For some satellites the NASA Prediction Bulletins are the only 
means available to the general user for obtaining orbital information. 
This paper provides the user with a computational interface between 
the information given in the NASA Prediction Bulletins and stanoard 
orbit determination programs. Such an interface is necessary to ob- 
tain accurate orbit predictions. The theoretical considerations and 
their computational verification in this interface morlelling are presented 
in detail. This analysis was performei in conjunction with satellite 
aided search and rescue position location experiments where accurate 
orbits of the Amateur Satellite Corporation (AMSAT) OSCAR -6 and 
OSCAR-7 spacecraft are a prerequisite. 


ORBIT MODKLLING FOR SATELLITES 
USING THE NASA PREDICTION BULLETINS 


I. INTRODUCTION 

Recent feasibility experiments with the Radio Amateur Satellites OSCAR-6 
and -7 indicate that a craicial factor in these experiments involving distress 
beacon location accuracy is the method of orbit acquisition and prediction utilized 
by the NASA Prediction Bulletins.* As part of a proposed satellite aided Search 
and Rescue mission (SAR), the beacon latitude and longitude is to be determined 
by using information from individual passes of Doppler frequency shift measure- 
ments in combination with satellite orbital data (Reference 1). The OSCAR orbital 
data is determined by use of a simplified generalized pei'turbations program 
(SGP), in conjunction with tracking data provided by the Space Defense Center 
network (Reference 2). The SGP utilizes a simplified first order Kozai analjdic 
orbit theory in order to calculate position and velocity coordinates of the OSCAR 
satellites. 

In this paper, the expression for the semi -major axis of the oibit is 
rewritteti co as to include the full first order short period perturbations. These 
effects are then introduced into the position and velocity magnitudes, thus re- 
setting the sidereal dlock, making the epoch vectors of thp :L;tate commensurate 
with those obtained from the standard Kozai and Brou theories. 

An incomplete description of the second order short period perturbation in 
the semi-major axis, produces no significant error in the Kozai nean 


*Scc “lixplanaiinn of The NASA Prediction Bulletins,” distributed by Operations Center Branch, Code 
512, NASA Goddard Space Fliglit Center, Grccnbelt, Md. 20771 (February 1974). 




mean-motion when the SGP orl>it constants arc chosen to proi uce a best fit to 
observations. However, the neglect of any terms of the first order, and often 
in the second order, gives a wrong instantaneous mean moticn for the start of 
a numerical integration, which then propagates as a secular like error in time. 
This requirement for the semi-major axis is peculiar to tha operation of using 
osculating elements to transfer the orbit computations to special perturbations, 
or even to another general perturbation. 

At present, there is a need at the National Aeronav.tics and Space Adminis- 
tration for a program which will reconcile the SGP results with the approach of 
the more sophisticated orbit theories /thereby permitting simultaneously accu- 
rate orbit determination and distress beacon location for considerably longer 
prediction arcs. 

Section II describes the SGP model and compares it to that of Brouwer and 
Koaai. Modifications to SGP are made in Section III. Sectio.ns IV and V describe 
the results and conclusions of calculations and comparisons of the modified SGP 
with the Space Defense Center SGP and the Brouwer theory. 

II. THE MODEL 

Brouwer, (Reference 3) and Korai, (Reference 4) have described the motion 
of an artificial Earth satellite without air drag. 

In the approach of Brouwer, the equations of motion of a small mass at- 
tracted by a spheroid are written as: 

X _ 3U d^y _ 3U d^ z _ ^ 

dt* Sx' dt^ By’ dt2 Bz’ 





with 


6 


U 



, ^*‘4 / 

(1-3 sin^/S) + h 


10 sin^ /3 + 


35 

3 



... . ( 1 ) 


Here, the equatorial plane of the spheroid is taken as the Cartesian xy plane; (i 
is the latitude, and if M is the mass of the spheroid and k the Gaussian constant, 
n - k^ M. The Delaunay variables are Introduced: 


G 

H 


a) , 


L(1 - 
G cos I , 


f 

g 

h 


mean anomaly, 
argument of the pericenter, 
longitude of ascending node. 


( 2 ) 


where a and e are the osculating semi-major axis, eccentricity and inclination 
of the orbital plane with the orbital plane respectively. With these, the equa- 
tions of motion become, 



F Is the Hamiltonian. 
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The problem is solved by considering a canonical transformation from the 
variables L, G, H, f , g, h to a new set L', G', H', f, g', h', and from the old to 
a new Hamiltonian F* (L', G', H', p , h'). This is accomplished with the aid 
of a determining function S(L', G', H', t, g, h), which is chosen in such a manner 
that f ' is not present in F*, while g' is permitted to appear. Consequently, L* 
and H' will be constants, and the system is reduced to one of one degree of 
freedom: 


dG d g* _ 3F* 

dt 3g' ’ dt 3G* 


( 4 ) 


After this system is solved, (' and h‘ are obtained by quadratuires from 


li! - _ iE! ItL - ^ (5 

dt " “ 3L' ’ dt ” ” 3H' ■ ' 

The determining function is obtained by a method used by von Zeipel (Ref- 
erence 5), in a qualitative study of the motions of minor planets. The final 
equations of motion are: 


a - a 


"|l + 72 (- 1 + r'3 ~ ~ r'^ (2 g' + 2 f ') j- 


e = e' + Sje + 


7-[r2[(-i + -7,-3^ 


,+ 3 (1 - 03) 


-T7-^^cos(2g* + 2f')j 
~ Tj' ( 1 - 03) 1^3 e* cos (2 g' + f ' ) + e' cos (2 g' +3 f ' )"j 


I 

i 

I 


f 


Ni 


I = I" + Sjl + 2 7:.0 (I - [3 cos (2 g' + 2 f') 1 - 3 e" cos (2 g' + f’) 


+ e" cos ( 2 g' + 3 f ' )j 


e = r 




c f 

Sin f 


3(1 - ^ ~ ^ + l^sin (2g' + f') 


/a"2 , a" l\ 

, ~jsin(2g' . 3f') 


} 


- ' + ^ 
^ - g %e 


» >2 |2 (- 1 + 3£»^)^^ 77 ^ + ^ + 1 j sin f ' 


( 6 ) 


3(1 - ^.2 “ r' ^ l^sin (2 g' + f') 


+ jsin (2 g' + 3 f')J j- 

+ ~ 7 j -^6 (-l + 56'^)(f'-f'+e"sinf') + (3-50*)j^3sin(2 g'^ 2 f ') 

+ 3 e" sin (2g* + f*) + c"sin(2g' + 3f')J^ 



\ 
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h " U' ~ ^ y 2 e" sin f') - 3 sin (2 g' + 2 f ') 

- 3 e* sin (2 g' + f * ) ~ e" sin ( 2 g' + 3 f ' )J. 

f, r' are to be computed from 

E' - e' sin E' = £' 


tan 


1 . / 1 + e”\'^ 1 

2^=1 2 ^ 


a* _ 1 + e* cos {' 

r' 


- «"2 


1 - e 


or 


~ sin f* = (1 “ sin E* 


cos f' - cos E' ~ 


~ ” 1 - e" cos E' 


,, Kj 

Here, f is the true anomaly, tj = (1 - e*’) , 0 = oos I', 7 j = 

a * 

and kj is related to the coefficient of the Earth's second zonal harmonic Jj , by 


the expression. 


^2 ~ ^2 »> • 


2 


m 


where is the Earth's equatorial radius. 


’ , / 





The quantities in brackets are of shcrt period, while a*, e", and I" are 
secular in nature. Secular and long period terms of the mean anomaly, argu- 
ment of perigee, and longitude of the node are contained in f ', g' and h', while 
long period variations in eccentricity and inclination are given by SjC and Sjl 
respectively. The periodic terms, both long and short, are developed tc 0(k^. 
The secular motions are obtained to 0 {kj^>. 

For the calculation of the coordinates at any time, the complete values of 
e and f should be used for the solution of Kepler's equation, 

E - e sin E = f (S) 

and subsequently, the radius vector r, and the true anomaly f which may then be 
used in the formulas: 

X = r [cos (g + f ) cos h - sin (g + f ) sin h cos I ] 
y = r [cos (g + f ) sin h + sin (g + f ) cos h cos I ] 
z - r sin (g + f ) sin I . (9) 

In a similar method by Kozai, perturbations of six orbital elements of a 
close Earth satellite moving in the gravitational field are obtained as functions 
of mean orbital elements and time. No assumptions are made about the order 
of magnitude of eccentricity and inclination. It is assumed however, that the 
density distribution of the Earth is symmetrical with respect to the axis of 
rotation. Also, the coefficient of the second harmonic of the potential is taken 
to be a small quantity of the first order, while those of the third and fourth 
harmonics are of the second order. In addition, the expression of the semi- 
major axis contains no long periodic terms. 


f '.f 

4 


The results include periodic perturbations of the first order, and secular 
perturbations up to the second order: 


a - a + da , 




e - e + de^ - de^ + dcj, 

i = T + di^ - di j + di , , 


- Wq + it + dc<J^ - cla!„ + doij , 


n = + fit + dn^ - + dfii , 


M ~ Mg + n t + dM 


" " "o{l+ f sinM^ }/r^ |. 


The secular perturbations of the first order are 


5 . , 

„ sin" 1 t 


il - fig - n t cos i 


M - Mg + n t , 


n - Ho + - On f 1 “ „ sin* i ) lA - e* 


I 
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where 00 ^^, and My are mean values at the epoch, that is, the initial values, 
from which periodic perturbations have been subtracted, n^ is the unperturbed 
mean motion, which is related to the unperturbed semi-major axis Og by 


n„2 Uy3 = G M. 


Kozai chose to adopt as a mean value of the semi-major ajds not a but 


A./ 


1 - -( 1 
I P" \ 


Sin' 


A 'XT 


SO that the following relation holds: 


■ GM<1 


p> ('■ 2 


( 12 ) 


The transformation from Kepler elements to Delaunay's canonical ones is: 


L - j/^i a , G - j/;ta ( 1 ~ e^) , H = j/ /xa (1 - e^) cos i , 


f = M, 


R - 


h = fl, 


(13) 


where jj. = G M, G is the gravitational constant and M is the mass of the Earth. 
Aj is taken to be of the first order of small quantities and is equal to 3 of 
Brouwer above. The quantity p = a(l - e^). Also, e and i are mean values 
with respect to M and &j, and a>g, fly, and My are also initial values from which 
periodic perturbations have been subtracted. 

The subscripted f quantities are long period perturbations, and are given 
by: 


6 
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^2 e sin^ i fl4-15 sin^i ^4 18-21 sin^i 


de, = dc* + — 

P ^ 4(4-5 sin^ i ) 


^ 18-21 sin^ i 1 

'a^ 7 I 


cos 2a; 


3 ^3 

+ ~ • — sin 1 sin Co, 

4 Aj a 


^2 sin 2i f 14 -15 sin^ i ^4 18-21 sin^ i 1 

."S , " I r cos 2co 

8 (4-5 sin^ t ^ ^2 ^ J 


3 ^3 


e cos 1 sin co , 


«2 e 

df), = dn - — — ^ 
’ p2 2(4- 


^ cos i 7-15 sin* i ^4 

• 5 sin^ i ) 6 ^2 


9-21 sin* i 


2 7 


5 sin* i ( 14 - IS sin* i ^4 18-21 sin* i 

+ ■< - — • sin 2co 

2(4-5 sin* i) 1 ® ^ 


3 ”3 cos i 

+ T e cos to , 

4 Ajp sin 1 


3 ^2 , . 

dto, - dto - ~ — sin* i sin 2-'o 

‘ ‘ 8 p2 


^2 1 fl4 - 15 sin* i 

_ * 

j* 4-5 sin* i I 


28-158 sin* i + 135 sin^ i 
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A 4 / I8-2I si 
A/ \ 28 


36 -210 sin^ i + 189 sit/ i 


^ sin^ i (13-15 sin^ i) / 14 - IS sin* i 18-21 sin^ i 


(4-5 sin^ i ) 


24 A| 28 j 


sin 2co 


3 ^3 sin^ i - cos* i 1 

+ ~ , — ^ ^ ~ cos Ol. 

4 A,p sin 1 e 


The barred quantities are those short periodic perturbations with respect 


to the mean anomaly, 


A 2 1 ~ 6 ^ - ■ ■■ 

de, = sin* i ~ cos 2v cos 2o>, 
* p2 6e 


^2 / 1 1 - e * \ 1 I 

do) = — ■< sin* i 1 r + cos 2v 1 + “ cos* i cos 2v ?• sin 2ai, 

‘ p* 1 \8 6e* ye 


1 ^2 


12 p2 


— sin 2i cos 2v cos 2o>, 


- - — — cos i cos 2v sin 2u>. 


/ 
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Here, 

cos jv = [ / : ) (1 + j - o^). 

\ 1 + V I ~ e 

where j = 1,2 The mean values of perturbations are not zero, except for 

the semi-major axis. 

The quantities it and oj for the node and argument of perigee account for 
secular changes through the second order: 



Li 


T 
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whore p ~ a (1 - e*), and! and e are mean values of the inclination and the 
eccentricity over all periods. The leading terms of these expressions aecount 
for the first order secular variations in H and u>. Both the long period and 
secular perturbations ai'e taken through second order. 

The short periodic terms contain terms only of the first order however. 
The results for the six elements are as follows: 


da, = “ [~^1 - ^ sin^ - (1 - j + sin^ i cos 2 (v + ^j)j. 




- „2s-3/2 


(1 - e^y 


~ (f) 2 (v + w)j 


sin* i ^ 


i ^ j 

ap 1 


cos 2 (v + oi) + e cos (v + 2o>) T c cos (3v + oj) 

o 


A, f / 5 \ 

dw, = — j [ \ ^ ~ ^ iy (v - M + e sin v) 

+ ^1 - ~ sin* I ^ ” 4 sin V + ~ sin 3v j- 

- ^ sin* * ^2 " M 


i ■ <’ 
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c 

+ — s>n 
16 


in^ i sin (\- “ 2c<j) - ^ 1 - sin^ sin 2 (v ^ w) 




i ) f sin (3v + 2co) 


+ sin^ i sin (4v + 2co) + sin* i sin (5v + 
8 lo 


di. = T “T sin 
‘ P* 


2i I cos 2 (v + w) + c cos (v + 2w) + “ cos (3v + 2 &j) | , 


jn = - — cos 

* p" 




M + e sin " 'J sin 2 (v + w) 


- sin (v + 2co) - — sin (3v + 2o)) 


‘ P* 



1 - I sin* i) I (^ ■ sin V + | sin 2v + ^ sin 3v| 


+ sin* i 


I A ^ ~ 16 ~ 


- ^ ^ jsin (v + 2w) - ” e sin (4v + 2u>) - ^ sin (5v r 2w)| j 


V 
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Here (v + cJ) is the true latitude, u. 

The SGP is a general perturbations program whose basic theory is that of 
Kozai. Its gravitational part is truncated to include only the effects of the first 
three zonal harmonics of the Earth's potential. In addition, the final equations 
of motion are carried through the first order accuracy in terms of the small 
parameter Jj, the coefficient of the Earth's second zonal harmonic. Included 
in the algorithm is a drag formulation which is based on the assumption that a 
mean element set will have its time variation given by the first few terms of a 
Taylor series. 
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The olassical mean element input set includes the inclination i^ , right 
ascension of the ascending node fig, eccentricity Cq, argument of perigee 
mean anomaly M^, and mean motion After updating the mean elements for 
first order seculai r.nd long period effects, the following quantities are calcu- 
lated, which involve the position, velocity, and true latitude respectively: 

r - a ( 1 - e cos E) 

• - . p 

r “ e sin E 




sir (E + oj) 


*yNSL " ^xNSL 


e sin E 


1 + 




and, 


a r 
< 
r 


cos u - “ cos (E + a?) - a 


xNSL 


e sin E ~\ 

®yNSU • 

1 + 


( 21 ) 


(E + cu) is obtained from Kepler's equation by an iteration procedure, a and e 
include atmospheric perturbations, 

Pl = a (1 - 




“xNSL 


C COS Ct>, 


SO 


II 



®yNSL “ ® 


2 Ja a (1 - e^) 


The subscript, SO, and L, denote secular and long period effects, and , the 
magnitude of the Earth's equatorial radius. 

The planetocentric distance, true latitude, ascending node, and inclination 
are then eorrected for short- period perturbations as follows: 


r = r + T !■) sin^ i- cos 2u 
4 ■* Pl ° 


u - u “ o Ja — r ^0 *■ 

O r\* 


^ ~ ^so + 4 J 2 ~ ^0 si" 2u 


i - in + 


sin in • cos cos 2u 


From these quantities, the osculating position and velocity can now be computed 
for a specified time. 

I 

The SGP is thus a first order analytical approach which includes perturba- 
tion effects through the Eartli's third zonal harmonic. Although the theory is 
not singulai* for any elliptic orbit, the position prediction is normally within 
±5 kilometers of the actual position between updated mean elements. It is this 
point to which we address ourselves in the next section, in order to make SGP 
consistent with both the Kozai and Brouv. or approach. • 


r 
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III. FIRST ORDER SHORT PERIOD FORM OF THE SEMI-MAJOR AXIS 
Adopting the form of the force function used by Vinti, (Reference 6), the 
last, and dominant term of the short periodic contribution to the semi-major 
axis (equations 17), can be written. 



For short period effects, we can use the relation given by Tisserand (Refer- 
ence 7), 




(1 - 


and consequently, equations (23) and (17) become. 


da 


S 


A 

2a (1 - 


Sill' 


cos 2u] . 


(24) 


Let us now consider the short period contribution to the planetocentric 
distance as given by equation (20), If as a first approximation, we assume the 
short period effects in the eccentricity to be small, due to the presence of the 
1/a* factor, then this parameter Is composed of the mean, secular, and long 
period effects, and let us designate this as e^^. In addition, the eccentricity can 



22 


be considered as constant over short intervals of time. We then have that 

da^ • (1 - ej) 

(1 + cos v) 

Over a period of revolution we have also that 


1 r ^ . 2- 

r 1 


4 

2n (1 + cos v) 2 77 

.Jo CCS v)J 


.j 

t 

■j 

1 

• (1 > e > 0). 

(26) 


Vl - 




Substituting this result and equation (24) into (25), we obtain 


f 

3 



.5 

= I ^ ^ ' 

cos 2u, 

(27) 



as the short period contribution of the semi-major axis to the planetoccntric 
distance of the satellite. This result differs from the first equation of (22), in 
that the leading coefficient here is 3/2 instead of 1/4. This discrepancy, and 
hence that of the corresponding semi-major, produces no significant change in 
the SGP mean mean-motion when the orbit constants are chosen to produce a 
best fit to observations. However, the factor of one quarter produces an Incorrect 
instantaneous mean motion for the start of a power series or numerical Inte- 
gration, and a corresponding discrepancy in the mean mean-motion results. 

Thus, equation (24) has the effect of resetting the sidereal clock for the satellite 
orbit. 


1 ' I ! 
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To further modify the SGP, the changes in the semi-major axis and 

planetocentric distance are then incorporated into the expressions for r and r v 
of equations (21), by the inclusion of the short period effects in the velocity. 

The velocity components are now calculated after the orbital parameters of 

equations (22) are determined. 

Following the approach of Kozai, and paralleling equations (10) and (11) of 
Section II, since the mean value of the semi-major axis is taken to be 


a = 



then the following relation holds: 


from which. 


n - 



Here, 



(28) 


(29) 


(30) 


where n^ is the unperturbed or mean mean-motion provided by the two-line 

element input. In our modified version of the SGP, we first calculate *>.„ from i 

I 

; i 

i .j 

. ..ti 
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equution (30) above, and then ntilize a and n from equations (28) and (29) together 
as starting values for the mean semi-major axis and mean motion, whereas the 
Space Defense Center SGP program utilizes a from equation (28) in conjunction 
with the given value of the mean mean-motion, n^. 

IV. RESULTS 

Table I is a comparison between the Space Defense Center SGP and 
Goddard (932) SGP with a Brouv. er-Navigational Analysis Program (NAP) 
system involving real data from the OSCAR-7 satellite. This data Is fitted with 
Space Defense Center SGP analytic theory and published in the NASA Prediction 
Bulletin number 87. The epoch of the data arc for bulletin 87 is, November 2, 

16 hours, 47 minutes, 51.989856 seconds. The fitted orbital elements for this 
bulletin are as follows: 

Inclination, i = 

Right Ascension of the Ascending Node, fi = 

Eccentricity, e = 

Argument of Perigee, a> = 

Mean Anomaly, M = 

Mean Motion, = 

'iLo Space Defense Center SGP, the Brouwer and Kozai methods, as well as 
the Goddard J?2 SGP programs employ basically the same fundamental theo- 
retical approach to the satellite motion problem, as opposed for example, to 
that of Vinti (Reference &}. While both SGP programs and Brouwer and Kozai, 
all treat the i)eriodic variations through fir.st order, only the latter two account 
for second order accuracy in the secular motions. This contrasts to the highly 


101.6817 degrees 
351.4182 degrees 
0.0012030 
269.8988 degrees 
90.0660 degrees 
12.53311719 ^^volutipns 


i 

I 
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accurate, non perturbation, separable Hamiltonian approach of Vinti which 
carries all terms to a minimum of third order accuracy. In addition, both the 
mean-motion and the short period form of the semi-major axis differ between 
Space Defense Center and Goddard version of the SGP. In the latter program, 
these quantities are expressed similar to those of Kozai, and in the case of the 
OSCAR 6 and 7 satellites, this difference in the semi-major axis, yields an 
initial discrepancy of approximately eight kilometers in position between the 
Space Defense and Goddard 932 SGP at the epoch. 

The primary differences between the Space Defense Center and Goddard 
932 SGP methods can essentially be attributed then to the respective definitions 
of the mean-motions, the corresponding mean and osculating semi-major axes, 
and the resulting short period contribution to the planetocentric distance of the 
satellite. 

Tables lA and IB show the comparison of Cartesian coordinates at the 
epoch of NASA Prediction Bulletin number 87. Differences in these coordinates 
between the Space Defense SGP, the Goddard 932 SGP, and the Brouwer analytic 

theory can be attributed to several causes. Tlie Space Defense program calcu- 
lations were fitted to the OSCAR 6 and 7 data. Since the two SGP programs are 
fundamentally similar to the theory of Brouwer, then there is sufficient reason 
to warrant a second order term In the short period perturbation of the semi- 
major major axis of the SGP version, as v/as done by Cohen and Lyddane 
(Reference 8) for the Brouwer theory. 

This ifi presently being undertaken by the authors. In addition, there is the 
added sophistication of the Brouwer and Kozai methods, that is not present in 
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either of the SGP calculations. In the former, the secular motions are developed 
to second order, while in SGP, these terras are carried to first order only. The 
comparison after three days is somewhat more striking, with the effect of the 
mean-motion clearly evident. These results arc listed in Tables IC and ID. 

Very similar results were obtained for the NASA Prediction Bulletins 
96 and 276, and these are shown in Tables II and III. Bulletins 87 and 96 are 
Specific epochs for OSCAR-7, and bulletin 276 relates to the OSCAR-6 
spacecraft. 
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Table lA 


Comparison of Cartesian Coordinates at the Epoch 
of NORAD Bulletin Number 87 



NORAD SPG 

Goddard 932 SGP 

Brouwer 

X = 

7743.758 km 

7748.970 knj 

7748.569 kmi 

y = 

-1168.462 km 

-1170.767 km 

-1169.135 km 

Z = 

-0.739 km 

6.515 km 

-0.979 km 

k = 

-0.213 km/s 

-0.216 km/s 

-0.213 km/s 

y = 

-1.428 Imn/s 

-1.427 km/s 

-1.428 km/s 

Z = 

6.987 lun/s 

6.984 km/s 

6.984 km/s 


N 



Table IB 


Ax « 
Ay = 


Differences of Cartesian Coordinates at the Epoch 
of NORAD Bulletin Number 87 


Brouwer-Is'ORAD SGP 


4.816 

km 


-0.673 

km 


-0.240 

km 


0.000 

km/s 


0.000 

km/s 


-0.003 

km/s 



Brcuwer-Goddard 932 SGP 


-0.401 km 


1.632 km 


-7.494 km 


.003 km/s 


-0.001 km/s 


0.000 km/s 













Table IC 


Comparison of Cartesian Coordinates of the Propagated State After 3 Days 
For NORAD Bulletin Number 87 



NORAD 

SGP 

Goddard 932 SGP 

NAP* 

X = 

-2404.339 

km 

-1589.378 

km 

-1586.169 

km 

y = 

-1284.918 

km 

-1406.800 

km 

-1408.339 

km 

z = 

7340.679 

km 

7533.466 

km 

7539.555 

km 

X = 

-6.753 

km/s 

-6.950 

km/s 

-6.946 

km/s 

y = 

1.097 

km/s 

0.965 

km/s 

0.964 

km/s 

z = 

-2J121 

km/s 

-1.288 

km/s 

-1.284 

km/s 


♦These values are the result of the numerical integration by the NAP program of tlie Brouwer epoch vectors. 
Tlicse results agree with the propagated values provided by both the Vinti ar i Brouwer analytic orbit 
theories. 
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Table ID 

Differences of Cartesian Coordinates of the Propagated State After 3 Days 
For NORAD Bulletin Number 87 


1 

NAP-NORAD SGP 

. J 

NAP-Godd.ard 932 SGP 

Ax = 

818.170 

Ian 

3.209 

km 

Ay = 

-123.421 

km 

-1.039 

km 

Az •■= 

198.876 

km 

6.089 

km 

Ax = 

-0.193 

km/s 

0,004 

km/s 

AS' = 

-0.133 

km/s 

-0.001 

km/s 

A z = 

0.737 

km/s 

0.004 

km/s 


\ 

} 

1 


& 

















Table IIB 


Differences of Cartesian Coordinates at the Epoch 
of NORAD Bulletin Number 96 

Brouwer-NORAD SGP Brouwer-Goddard 932 SGP 


-1.112 

km 

-7.773 

km 

-1.326 

km 

-4.807 

!un 

0.222 

km 

-12.644 

km 

0.001 

km/s 

0.005 

km/s 

0.000 

km/s 

0.004 

km/s 

0.004 

km/s 

9.008 

km/s 
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Table IIC 

Compai'ison of Cartesian Coordinates of the Propagated State After 1 Day 
For NORAD Bulletin Nu Tiber 96 



NORAD 

SPG 

Goddard 932 SGP 

NAP* 

X = 

64.828 

km 

276.813 

km 

266.782 

km 

y = 

-2294.326 

km 

-2066.497 

km 

-2078.811 

km 

Z = 

7470.485 

km 

7530.580 

km 

7532.645 

km 

X = 

-4.783 

km/s 

-4.776 

km/s 

-4.773 

km/s 

y = 

-5.087 

km/s 

-5.166 

kra/s 

-5.161 

km/s 

i. = 

-1.526 

km/s 

-1.250 

km/s 

-1.256 

km/s 


•These values are the result of the numerical integration by the NAP program of the Brou ver epoch vectors. 
These results agree with the propagated values provided by both the Vinti and Brouwer analytic orbit 
theories. 


Table HD 


Differences of Cartesian Coordinates of the Propagated State After 1 Day 
For NORAD Bulletin Number 96 



NAP-NORAD SGP 

NAP-Goddard 932 SGP 

Ax = 

201.954 km 

-10.031 km 

Ay = 

215.515 km 

-12.314 km 

Az = 

62.160 km 

+2.065 km 

Ax = 

0.010 km/s 

0.003 km/s 

Ay = 

-0.074 km/s 

0.005 km/s 

Az = 

0.270 km/s 

-0.006 km/s 














X 
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Table niB 


Differences of Cartesian Coordinates at the Epoch 
of NORAD Bulletin N imber 276 



— 

Brouwer-NORAD SGP 

Brouwer-Goddard 932 SGP 

. 

Ax = 

-4.979 

km 

-12.149 

km 

Ay = 

-4.078 

km 

-2.978 

km 

Az = 

13.440 

km 

-0.213 

km 

Ax = 

-0.005 

km/s 

+0.001 

km/s 

t> 

'<• 

II 

-0.003 

km/s 

-0,003 

km/s 

Az = 

0.007 

km/s 

0.011 

km/s 


/, . 
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Table niC 

Comparison of Cartesian Coordinates of the Propagated State After 6 Days 
For NORAD Bulletin Number 276 



NORAD 

SGP 

Goddard 932 SGP 

NAP* 

X = 

6828.173 

km 

7277.146 

km 

7267,588 

km 

y = 

1644.624 

km 

2118.325 

km 

2113.978 

km 

Z = 

3464.997 

km 

2010.906 

km 

2013.454 

km 

X = 

-2.625 

km/s 

-1.306 

km/s 

-1.311 

km/s 

y = 

-2.274 

km/s 

-1.921 

km/s 

-1.924 

km/s 

Z = 

6.233 

km/s 

6.739 

km/s 

6.746 

km/s 


•These values are the result of the numerical integration by the NAP program of the Brouwer epi.ch vectors. 
Tliesc results agree with the propagated values provided by both the Vinti and Brouwer analytic orbit 
theories. 
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Table HID 

Differences of Cartesian Coordinates of the Propagated State After 6 Days 
For NORAD Bulletin Number 276 



NAP-NORAD SGP 

NAP-Goddard 932 SGP 

Ax = 

439.415 km 

-9.558 km 

Ay = 

469.354 km 

-4.347 km 

Az -- 

-1451.543 km 

2.548 km 

Ax = 

1.314 km/s 

-0.005 km/s 

Ay = 

0.350 km/s 

-0.003 km/s 

Az = 

0.513 km/s 

0.007 km/s 



) 

k 


I 

I 


I 





i 


I 
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V. CONCLUSIONS 

The inclusion of the full first order short period effects in the semi-major 
axis and corresponding planetocentric distance in the SGP analjlic theory, gives 
a good instantaneous mean motion for the start of the Navigational Analj'sis 
numerical integration program. It has also been shown by Lyddane and Cohen 
(Reference 8), that in the case of a program such as Brouwer, the operation 
of using osculating elements to transfer the orbit computations to special per- 
turbations, or even to another general perturbation, has the additional require- 
ment for the semi-major axis to be of second order in short period terms. 
Hence, such a transformation would be of great value in eliminating any secular 
discrepancy. 

Finally, the mean value of mean motion, as well as the semi-major axis 
are redefined in accordance with the relationship given by Kozai. The resulting 
modified version of the SGP is found to agree much more closely both at epoch 
and in the prediction mode, with the Brouwer Anal 5 ’tic Theory and the Naviga- 
tional Analysis (numerical integration) Program. 

VI. ACKNOWLEDGEMENTS 

Mr. John S. Watson of the Goddard Space Flight Center carried out some 
of the numerical computations. Mr. R. A. Gordon of the Goddard Space Flight 
Center took part in several useful discussions about the work. 



Table HIB 


Differences of Cartesian Coordinates at the Epoch 
of NORAD Bulletin N imber 276 



Brouwer-NORAD SGP 

Brouwer-Goddard 932 SGP 

Ax = 

-4.979 km 

-12.149 km 

Lsy - 

-4.078 km 

-2.978 km 

Az = 

13.440 km 

-0.213 km 

Ax = 

-0.005 km/s 

+0.001 km/s 

Ay = 

-0.003 km/s 

-0,003 km/s 

Az = 

0.007 km/s 

0.011 km/s 





) 
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Table niC 

Comparison of Cartesian Coordinates of the Propagated State After 6 Days 
For NORAD Bulletin Number 276 



NORAD 

SGP 

Goddard 932 SGP 

NAP* 

L ... 

X = 

6828.173 

km 

7277.146 

km 

7267.588 

km 

y = 

1644.624 

km 

2118.325 

km 

2113.978 

km 

z = 

3464.997 

km 

2010.906 

km 

2013.454 

km 

X = 

-2.625 

km/s 

-1.306 

km/s 

-1.311 

km/s 

y = 

-2.274 

km/s 

-1.921 

km/s 

-1.924 

km/s 

Z = 

6.233 

km/s 

6.739 

km/s 

6.746 

km/s 


•These values are the result of the numerical integration by the NAP program of the Brouwer cpi>ch vectors. 
Tiiesc results agree with the propagated values provided by both the Vinti and Brouwer analytic orbit 
theories. 
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